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Abstract 

We present a study of the chaotic behavior of the bouncing ball 
billiard. The work is realised on the purpose of finding at least certain 
causes of separation of the neighbouring trajectories. Having in view 
the geometrical construction of the system, we report a clear origin 
of chaoticity of the bouncing ball billiard. By this we claim that in 
case when the floor is made of arc of circles - in a certain interval of 
frequencies - a lower bound for the maximal Ljapunov can be evaluated 
by semianalical techniques. 

1 Introduction 

Deterministic features of transport has been studied in different problems 
[1] . These works has also shown that transport may be related to the chaotic 
aspects of the dynamics [2] . 

The idea of bouncing ball was studied in different problems where ana- 
lytical approximations [3] and comprehensive numerical works can also be 
found jl]. The analytical approximations have been shown the possible evi- 
dence of bifurcations while the numerical works presented chaotic regimes of 
the bouncing ball system. 

The bouncing ball billiard as a spatial extension of the one dimensional 
bouncing ball problem has been introduced in [5]. This first work enhances 
the irregular diffusivity of the system which is similar to certain models 
of transport [6]. The following work has outlined the spiral modes in the 
phase space of this problem [7] and also pointed out its relevance on granular 
matter [8]. Idealised versions where the bounces are performed without loss 
of energy, i.e. the restitution coefficient is one, and there is no oscillation of 
the floor may be found in [9]. 



The chaoticity of the sawtooth type of the bouncing ball billiard has been 
studied in [TO], with considerable theoretical background [UJ. Considering 
the problem as a gravitational billiard, aspects on chaotic features are also 
approached by numerical methods in [12J. 

The present work focuses on the geometrical origin of chaoticity of the 
bouncing ball billiard where it is investigated the impact of the curvature of 
the arcs of circles on the maximal Ljapunov exponent. The derivation shows, 
that in case of resonance one may give semianalitical estimates on chaotic 
behavior. 

The study on manifolds for multi-dimensional billiards related to geomet- 
ric properties is presented in [T5] . 

From the practical point of view the reaction of CO with O2 on Pt sur- 
face, which is under thermal excitation, the molecules CO performs diffusive 
motion on the surface before the reaction would occur |14j . 

Quasi-deterministic aspects on diffusion may occur in the behaviour of of 
certain species where in the process of food searching one can find a random- 
ness, but there is also a kind of determinism because the animals may have 
certain remembrances on the places where they found food in the past [15J. 

The article is organised as follows. In Section [2] we shortly describe the 
bouncing ball billiard system. The Section [3] makes a presentation of that 
frequency region where the semianalical approaches to some extent are possi- 
ble. Section H] shows an evaluation which has a semi-empirical and analytical 
background giving a lower bound for the maximal Ljapunov exponent. Fi- 
nally, Section [5] discusses the similarities and differences of the analytical and 
real value of the maximal Ljapunov exponent. 

2 The bouncing ball billiard 

At this point we make a review of the most important features of the bounc- 
ing ball problem. The system studied is a point particle, which bounces on 
a floor realised of arc of circles |16j . The floor is oscillating with a frequency 
/ corresponding to a circular frequency u = 2nf. The system is presented 
below. The bouncing ball billiard from the point of view of the diffusion 
was presented in a comprehensive way in [5j [7] . The main conclusions are 
that the system possesses irregular diffusion, and the principal maximas for 
the diffusion occurs at the resonances. These resonances are at the frequen- 
cies, where the time of the flight becomes equal or multiple of the period of 



2 




Figure 1: The bouncing ball billiard. A point particle realises bounces on a 
vibrating floor consisting of arcs of circles. 



vibration applied. 

The bouncing ball billiard that we study in this paper, with the floor 
formed by circular scatterers, is depicted in Fig. [TJ 

The equations of motion of this system are presented as follows: The 
particle performs a free flight between two collisions in the gravitational field 
g || y. Consequently, its coordinates (x~ +1 , y^+i) and velocities (vZ n+ i,v~ n+1 ) 
at time t n+ \ immediately before the (n + l) th collision and its coordinates 
y+) and velocities (v£ n , u+ n ) at time t n immediately after the n th collision 
are related by the following equations 

X n+1 = X n + V tn(tri+1 — t n ) (1) 

Vn+l = Vt + U yn(Wl - *n) ~ #(Wl ~ t n ) 2 /2, (2) 

V x,n+1 = V xn (3) 

V y,n+1 = V+ n ~ g{t n+1 -t n ) . (4) 

At the collisions the change of the velocities is given by 

V ln - V ci±n = k (V ci±n - VlJ (5) 

P (v7 n - Vd\\ n ) , (6) 



V\ „ — v 



n u ci\\n 

where v C i is the velocity of the corrugated floor. We distinguish between 
the two different velocity components relative to the normal vector at the 
surface of the scatterers, where the scatterers are represented by the arcs of 
the circles forming the floor. v±, vu and v C i±, v c q is the normal and tangential 
components of the particle's, respectively the floor's velocity with respect 
to the surface at the scattering point. Correspondingly, we introduce two 
different restitution coefficients k and (3 that are perpendicular, respectively 
tangential to the normal. 

As in case of the vertically bouncing ball problem we assume that the floor 
oscillates sinusoidally, y ci = —A sin(a;£), where A and uj are the amplitude 
respectively the frequency of the vibration, see FigJTJ 
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The radius of circles are i?=15mm and the restitution coefficients k = 0.7, 
respective = 0.99. It is important that the slope on the arcs of the circles 
is very shallow. The distance between two arc of circles is d = 2mm. By this 
terms proportional with d 2 /R 2 or less are considered terms with second - or 
higher - order. 

2.1 Considerations on the chaoticity of the bouncing 
ball 

The chaos of the bouncing ball which spatially is a one dimensional system 
and which is performed on the vertical direction has been discussed in [I]. 

If we consider Ax the initial displacement and Ax n the displacement 
after n bounces between neighboring trajectories for the one dimensional 
vertically bouncing ball, then the Ljapunov exponent may be evaluated: A ~ 
(1/tn) log(Ax n /Axo). The main problem is that at certain frequencies the 
ball may be stucked on the surface. Because the neigbouring trajectory 
usually is also stucked, i.e. Ax n becomes zero, consequently in such cases 
the Ljapunov exponent A is simply undefined. 

Of coarse one may discuss on chaotic regimes between two stucks, but 
this is the reason why in general considerations about the chaoticity of the 
vertically bouncing ball should be done in a very careful way. 

In spite of the fact that in two dimensions it is almost impossible to have 
neighboring orbits stucked at the same place, the following study will try to 
avoid that frequency regions where stucking orbits might be possible. 

3 General considerations 

We discuss the chaoticity of the bouncing ball billiard at the frequency re- 
gion where the 1/1 resonance holds. The approximation we try to make is 
semi-empirical and the considerations are presented below. The first obser- 
vation we make is, that the time of the flight between subsequent collision 
at 1/1 resonance is approximately the same. We note this time with tfi y 
and corresponds to that time while the particle makes one bounce and it is 
close to the time while the floor makes one complete oscillation. Subsequent 
values of tfi y at 1/1 resonance are shown below 
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30 


0.0185 


31 


0.0190 


32 


0.0184 


33 



One can see that these values do not differ too much. They cannot be 
the same, for instance because the surface is not flat, but they are close to 
each other and around a specific value. 

The other observation is, that the first resonance manifests so that the 
elongation almost reaches its maximum A, and the particle meets the floor 
for almost all cases very close to this height A - see fig. EJ 

4 Geometrical origin of chaoticity 

We consider one component of the velocity, namely the horizontal one v x . 
We are interested in change in difference of the horizontal component of 
the velocities of two neighboring trajectories. By this we try to make an 
estimate of a lower bound of the Ljapunov exponent, which manifests on 
the v x direction of the phase space, which would also give a picture on the 
horizontal chaoticity of this problem. 

The figure below fig. |2] shows two trajectories, starting from the same 
place, with slightly different velocity vectors. The picture is at 1/1 resonance, 
where after a certain transient the bounces are made a little bit above the 
vertical coordinate A = 0.1 mm which denotes the amplitude. 

We assume that at the starting point there is a difference in angles but 
not in the magnitude of the velocities vq- This initial deflection in angles we 
denote by S. The corresponding difference in initial velocities we denote by 
Avx ini- Because of this initial deflection in angle there will be at final arrival 
a change in horizontal coordinates Ax. 

Due to the convex curvature of the arcs characterized by the radius R, 
the displacement Ax at the arrival on this curvature will cause further de- 
flection in angles after one collision which we denote by 5'. This cause a final 
difference in the horizontal component of velocities after the first bounce 
Av xJin . 

The rate of exponential separation of the trajectories after such a bounce 
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Figure 2: Typical trajectory of the 1/1 resonance. While the particle arrives 
at a height close to 0.6, horizontally in general a length around 0.15 is made. 
The first conclusion is, that the particle in most of the cases will arrive in 
a steep angle, so the incident angle relative to the normal to the surface is 
small. The other conclusion is, that quite a number of such bounces occurs 
while the particle arrives from one arc to another one. The second trajectory 
is a neighboring one, caused to a deflection of angles of initial velocities at 
the starting point. 



due to the finite radius R we denote with \ VXy R 

1 , \Av xJin \ 



XvxR = —\n " (7) 

tfly \A Vx : 



J x,vni | 



At the launch of the trajectory we consider the magnitude of the velocity 
v , the angles relative to the vertical are and a^ 2 \ the complementer 
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angles are and where as it was mentioned above cr 2 ) = + S 



I Au, 



\v cos a)p—VQ cos \ 



\v Q sin a 



(i). 



-f o sin a 



(2)1 



|f o 5 cos | +h.o.t 



where the higher order terms means terms which are proportional with at 
least the second power of 5. Correspondingly if 5' is the angle between the 
directions of trajectories after the first bounce 



I Aw 



x,fin | 



\vq 5' cos I + h.o.t 



(9) 



Because there is a free flight in gravitational field and the arcs are with shal- 



low slope practically a 



(i)' 



a 



(i) 



very precisely their difference is a second 



order term. We mention, that ar« is also small as one can see on the fig. [2J 
so its product with 5 is also considered a value with second order. 

As a result we get for the value X VX} r, which has a definite contribution 
due to the finite value of 



A 



vx,R 



In 



l fiy 



+ h.o.t. 



(10) 



where one can see that the ratio between the deflection of the angles after 
the bounce and before the bounce counts. 
The bounce is presented on Fig. [3j 

Based on Fig. [3] we can conclude that the deflection between the two 
reflected trajectories (5') one hand is due to the initial deflection 5. For the 
reflected trajectories there is a contribution due to the curvature. If the point 
of incidence of the first trajectory is at 9 then the point of incidence of the 
second one is at 9 + d9. The incident angles relative to the normal differs by 
d6 and in addition the reflected angles - on the Fig. |3] - j[ and 72 has also a 
difference dQ. Consequently the term d9 has to be counted twice. 



5' = 5 + 2d9 
Inserting this relation in eq.f llOp one gets 



t 



In 



fly 



S + 2d9 



1 , 

+ h.o.t. = In 

tfiy 







+ h.o.t. 



(11) 



(12) 



1 Even if a term proportional with a or sin a would be kept, at the end where the 
average value is calculated for X vx or it drops out or it is proved to be of higher order. 
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Figure 3: Illustration of the bounce of neighboring trajectories on an arc 
of circle. The figure enhances the displacement of trajectories and their 
velocities after the bounce due to the geometry of the floor. 

At the arriving point on the surface the initial difference 5 will cause a 
displacement Ax. This means that on the arc of the circle the trajectories 
will arrive at a difference d6 ~ \Ax\/R as one can see on fig. [3J Now follows 
the evaluation of Ax and correspondingly of the d9. 

The length Ax is due to the difference in angles of the velocities, at the 
starting point. The first one is launched at an angle the other one with 
an angle = a)P — 5. At the end one of the particles arrives at x 2 , the 
other one at x\ 

^_ 2 fK«-co 8 « (13) 

where is the angle of the velocity made with the horizontal direction of 
the first trajectory at the starting point. During the evaluation we make the 
approximation, that sin 5 is approximately 5 and cos 5 ~ 1 or the differences 
are at least second order in 5, and are included in the higher order terms. 

v 2 

Ax = x l -x 2 = 2-^5[cos 2 a™ - sin 2 a£>] + h.o.t. (14) 

At this point one can see that cos 2 a^ 1 ' = sin 2 - here being the 
incident angle relative to the vertical - can be neglected, consequently d9 
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yields the following value 



d 0^\M^ 2 S^-(sm 2 a^) (15) 
R gR 

If we take into account that the time for the flight is tji y = 2t> sin aj^'/g then 

2^4^(sin 2 a«)^^£ (i 6 ) 
5 gR R 

By this we get for the value \ VX) r in leading order 

\ m , R = ln(l + + h -°- L = 9 ~^ + h.o.t (17) 

tfly K K 

where the logarithm has been expanded, and terms proportional with 1/R 2 
have been also considered as being of higher order. 

The average of X VXt R means averaging the expression above. By this the 
higher order terms vanishes or becomes smaller so they remain of higher 
order. Consequently it reduces to the average of the time of the flight tfi y . 
Its average is given by the period of oscillations T resulting for X VXt R in 
leading order 

A ra , R - f - ^ (18) 
R J 

Because \ vx ,r is a manifestation of the separation of the neighboring trajec- 
tories in the v x direction due to the geometry, consequently this value can be 
considered a lower bound for the maximal Ljapunov exponent. 



4.1 The case of two periodic orbits 

As it is pointed out in the work |3j - with increasing the frequency - bifurca- 
tion of the resonant trajectory may be possible. This means that the time of 
flight consists of a shorter and a longer time alternating one after the other 
which we denote by and but they still do not differ too much from 
each other. @ 

2 In general in the case of the bouncing ball billiard to have a considerable difference 
between and tfi y ,i even it is not possible, or because the dynamics enters in further 

bifurcations, or simply it crashes to a scenario with lots of sticking orbits. 
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In such case the approximation that have been presented previously are 
still valid and one gets after two consequent flights - one is shorter, one is 
longer - 

A ra B = l' 1 



5 5' 



tfivl + tfiv,2 

The argument of the logarithm can be written as 



+ h.o.t. 



A 



1 



vx,R 



t 



fly,l + tfly,2 



In 



1 + 



R 



1 + 



R 



+ h.o.t. 



After the expansion to the first order we have 



A 



vx,R 



9_ 
R 



j.2 _i_ +2 



t 



flyA 



+ t 



+ h.o.t. 



fiy,2 , 



In the numerator of the second fraction the decompositions are made 



t 



fiy,m 



tfly,l + tfly t2 ± t 



Finally we arrive to the relation 

9 tfly,! + tflyj 



A 



vx,R 



R 



+ 



9 {tfi 



t 



fly,2) 



h.o.t. 



(19) 



(20) 



(21) 



(22) 



(23) 



R2\tfly,l + l fly,V 

This expression can be averaged and the term proportional with — 
tfiy^) 2 is still too small and is considered of second order. The average of the 
last expression yields the value 

A„, H .|PI±iy (24) 

where T\ and T 2 represents the average values of the respective t/zj/,2- 
Because even in the case of the two periodic orbit (7\ + T 2 )/2 equals an 
average time flight T' which is the inverse of that frequency / where the 
dynamics is already bifurcated. This is in fact an interval of frequencies, so 
the latter formula for two periodic orbits still shows a strong analogy with 
eq. ffl8l) and it may be written 



Kx,R — — 



9_ 

Rf 



(25) 



This latter formula shows that the relation ( fl8|) may be valid for the full 1/1 
resonance and for bifurcated trajectories not too far from it. 
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5 The real value of the maximal Ljapunov 
exponent 

In this section we discusses connections of the geometry with the chaoticity. 
The relation ( TT5]) in case of the 1/1 resonance at a given gravitational field g 
has a characteristic time of flight T. Because of this reason, one may consider 
that one of the most relevant dependence of the maximal Ljapunov exponent 
from practical purposes is in terms of the radius. So in this section we restrict 
ourselves to the case of the frequency / = 53.2 Hz when the first resonance is 
fully developed. B The following figure illustrates the numerical values of the 
maximal Ljapunov exponent in terms of inverse of the radius. Consequently 
we start at a curvature 1/24 mm -1 and we end at a curvature 1/16 mm -1 . 
The endpoints are also motivated by the works [5j[7]. 

60 | 1 1 1 1 1 



50 - 




40 - 
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T3 

| 30- 
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10 - 

o I ' ' ' 1 
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K 

Figure 4: Presentation of the maximal Ljapunov exponent as a function of 
the inverse of the radius, i.e. the curvature K, at the frequency / = 53.2 Hz. 
The continuous straight line presents the lower bound estimation f fT8|) . 

As one can see the maximal Ljapunov exponent increases as the curvature 

3 The interval where we expect a validity of the relation (|18|) is a frequency region 
/ = 53.2 ± 5 Hz. The validity sometimes may be even wider, however there is a slight 
dependence of the endpoints of the interval on R, consequently it has been chosen a 
common domain in which the resonance or eventually the bifurcated scenario holds for 
different radii we discuss below. 
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- the inverse of the radius - increases. This is normal because a smaller radius 
implies a stronger separation of neighboring trajectories. From quantitative 
point of view both the numerical and the analytical evaluations has similar 
behavior - as it is shown on the fig. HI 

Regarding the quantitative aspects the evaluation ffTB"]) is considerably 
below of the real value. Partly because the analytical evaluation is just 
a projection. On the other hand the evaluation wants to detect only the 
geometrical effects on the chaoticity. Of coarse considerable other effects 
may have an important role on the chaoticity, but the study of such aspects 
will need further work. 

We thank Rainer Klages for comments and encouragement. 
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